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Abstract 

The impact parameter and rapidity dependence of the Cronin effect for massless pions in d + Au 
reactions at \^s NN = 200 GeV at RHIC is computed in the framework of pQCD multiple elastic 
scattering on a nuclear target. We introduce a formalism to incorporate initial state energy loss in 
perturbative calculations and take into account the elastic energy loss in addition to the transverse 
momentum broadening of partons. We argue that the centrality dependence of the Cronin effect 
can distinguish between different hadron production scenarios at RHIC. Its magnitude and rapidity 
dependence are shown to carry important experimental information about the properties of cold 
nuclear matter up to the moderate- and large- x antishadowing/EMC regions. 
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I. INTRODUCTION 

The discovery of jet quenching at RHIC and its derivative signatures - the high- 

Pt azimuthal asymmetry [3] and the disappearance of back-to-back correlations ^| - are 
among the most direct evidence for the creation of dense QCD matter in nucleus-nucleus 
collision, Recent theoretic*. wor, gS has provided snccessfnl q ualitative and q ua„ tit a t ive 
explanation of the experimental data via the techniques of perturbative QCD with strong 
partonic final state interactions leading to non-Abelian energy loss and related modifications 
to the fragmentation functions . The important question of the global p^ and 

y/s systematics of the nuclear modification factors Rba{pt) m P + A and A + A reactions 
has also been addressed j^j]. 

Alternative mechanisms related to hadronic rescattering/absorption, initial state wave- 
function effects, and coherent nucleon scattering have been suggested [Ll| as possible ex- 
planations of the observed high-p^ hadron deficit. This letter extends the discussion of 
in light of the upcoming d + Au data and demonstrates how the centrality dependence of 
the nuclear modification factor can distinguish between contrasting physical pictures. We 
emphasize the need for comparison between data and theory at moderate and high pt since 
vastly different models give comparable results at the level the inclusive dN ch /dy. 

The importance of the d + Au measurements at RHIC, however, reaches far beyond test- 
ing the validity of the perturbative baseline for jet quenching calculations. In the transverse 
momentum 2 < px < 10 GeV and rapidity —3 < y < +3 intervals considered here the mea- 
sured hadrons retain a memory of the nuclear modifications to the parton distribution 
functions (PDFs) in the 10~ 3 < x < 1 range. We find that in the small- x region of the 
nucleus the multiple elastic scatterings of the incoming partons completely dominate over 
nuclear shadowing. For large nuclear x-values the EMC effect and the initial state energy 
loss may lead to moderate (20-30%) suppression relative to the binary collision scaled p + p 
result over a wide px range. The detailed transverse momentum dependence of the nuclear 
modification ratio may provide valuable complementary information about the gluon anti- 
shadowing/EMC region, which is not well constrained by data. In fact, the x-dependence 
in the gluon shadowing function S 9 /a(x, Q 2 ) may be significantly weaker than currently 
anticipated. 
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II. SCALING OF HADRON SPECTRA IN d + Au REACTIONS 

The Cronin effect is defined as the deviation of the total invariant hadron cross 
section in p + A reactions from the binary collision scaled p + p result. For the more general 
case of B + A reactions the scaling factor is given by B ■ A. The original parametrization 
of dynamical nuclear effects via the Cronin power a{pr) relates to the nuclear modification 
ratio R B a(pt) as follows: R p a{pt) = A a ^~ l . The same relation (without additional factors 
of 2 since the partons from the nucleus do not scatter significantly on the deuteron) can also 
be used to extract a total inclusive cedAuipr) an d a centrality dependent a dAu (PT,b) from 



RdAu(P 



da dAu ,Ari ■ A A „ da pp 



I — ~ j J-? ' total invariant cross section, d + Au 



T) 



dyd 2 pT dyd 2 pT 

dN dAu (b) ,T dA Jb) da pp , , 7 , 

— — / : — , about impact parameter o, a + Au 

dyd z px dyd 2 px 

In Eq.flTJ) T d Au(p) = J d 2 r T d (r)T^(b — r) in terms of the nuclear thickness functions 
Ta(P) = J dz p^(b,z) and the scaling factors are computed according to the Glauber 
model Q]. It is important to note that both enhancement (pr ^ 1.5 — 2 GeV) and sup- 



pression (pt < 1 — 1.5 GeV) are an integral part of the Cronin effect. Experimentally 14], 
at pt = 0.77 GeV a{px) is below one, leading to a factor ~ 2 suppression. However, the 
deviation of a(pr) from unity is reduced with increasing yfs and decreasing xt = 2pr/y/s - 
a trend exactly opposite to the expectation from strong shadowing and/or gluon saturation. 
The pt positions of the Cronin peak R p a max and intercept R p a = 1 are also stable, i.e. they 
are strongly x^-dependent. 

The lowest order pQCD differential cross section for inclusive A + B — > h + X production 
that enters Eq. (JTJ) is given by j^, [3] 



1 da- 



dAu 



A d ■ A Au dyd 2 p T 

1 dN dAu (b) 
T d Au{b) dyd 2 p T 



K • dx a dx b / d 2 \s. a d 2 k b ® g(k a )g(k b ) ® S d (x a , Q 2 a )S Au (x b , Q\ 

abed J J 



® fa/d(Xa, Ql)f b/ Au(Xb, Ql) » ® D ^ M ^ , (2) 

dt 1TZ C 

where we have used the same sets of lowest order (LO) PDFs, fragmentation functions 
(FFs), and shadowing functions as in Q. In Eq. (J2J) the renormalization, factorization, 
and fragmentation scales are chosen to be \i = Q a ^ b = Q c = Pt- 



The the fcr-broadening function in (J2J) approximates generalized parton distributions: 

1 



fa/p (k T ; x a ,Q a ) ~ g(k T ) ® 

fa/p{%a, Qa) 



For illustration we include a discussion of parton broadening relative to the axis of propaga- 
tion in the framework of the leading double log approximation (LDLA). The k^ probability 
distribution resulting from vacuum radiation is given by 18] 



(3) 



1 da 



o" dk\ 



-2 
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2n ky Q 2 
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(4) 



with mean broadening (k|) pp = (7Ra s /7r(l + 0(a s ))Q 2 . Cr = {Cf,Ca} is the SU(3) 
Casimir in the corresponding fundamental or adjoint representation. Despite the strong mo- 
mentum ordering approximation , which gives an unrealistic null probability for retaining 
the original jet direction, Eq.(@]) qualitatively describes many the jet physics features. For 
example, with ~9 m = a/ (k^pp/Q 2 ~ a/ CrOs/it being the estimate for the half-width of the 
jet cone a slow logarithmic narrowing with Et is predicted. Similarly, to lowest order the 
cone of a gluon initiated jet is found to be 50% wider (a/C^/CV) than the cone of a quark 
jet. 
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(k 2 T ) pp [GeV] 
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1.000 
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1.082 



TABLE I: The fractional contribution of valance quark hard scattering to hadron production at 
y = 0, Q = pt — 2.5 GeV versus a/s. The vacuum radiation induced parton broadening computed 
from Eq.© and its variation A(y / s) relative to the a/s = 27.4 GeV reference are also shown. 

In this letter we go beyond the A/s-independent description of the radiative and elastic 
broadening and take into account the initial parton subprocesses for produced hadrons. An 
increased contribution of hard gluon scattering results in an enhanced 



04)pp = (CfF^VS, Q) + CUF*(V5, Qj) c 2 



(5) 



In Eq.<$ F q h (yfi,Q),F^,Q) 



Fq{y/~s, Q) are the fractions of hadrons coming from 



the hard scattering of valance quarks and gluons/sea quarks. With sea quarks we associate 
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the Casimir of their gluon parent (g — > qq). For hadrons with pt = Q — 2 — 3 GeV where the 
fc^-smearing is reflected, the fraction coming from valence quarks is computed from collinear 
factorized pQCD and given in Table|TJ The corresponding (k^) pp is found form Eq.© where 
C 2 = 0.646 GeV 2 has been fixed through comparison to the y/s = 27 A GeV data |l4] on 
i(7r+ + 7T - ) production in p + p. A( v /i) = (k^( v /i)) pp /(k^(27.4 GeV)) pp in Table □ directly 
reflects the variation with yfs of the effective mean color charge of ~few GeV projectile 
partons relative to the \/s = 27 A GeV case. 

We have checked that the energy and transverse momentum dependence of inclusive 
hadron production cross sections is well described with the computed values of (kf.) pp which 
exhibit a small < 10% variation about (kf,) pp = 1.8 GeV 2 used in j^j]. 

III. IMPLEMENTATION OF INITIAL STATE ENERGY LOSS 

While there are indications that the effect of initial state gluon bremsstrahlung may be 
small the elastic energy loss of propagating hard probes in cold nuclear matter has 
not been considered before. Discussion of QCD scattering in nuclear targets can be found 
in j^. In both frameworks of the Qiu-Sterman approach and the Gyulassy-Levai-Vitev 
reaction operator approach the longitudinal correction that corresponds to the transverse 
momentum broadening can be expressed through the momentum shift operators e 
The initial parton flux d 3 N l (k\\, kjj flux is propagated to the hard collision vertex according 
to bJl: 



rf 3 V/(fc||,k ± ) = yX 

dk\\d 2 ~k\ ^ n\ 



00 ,.n r n 



n=0 " i=l 



1 d(T ^L ( e -a^x® e -*Vi% 



Tel d 2 q i± 



d 3 N*(k lb k ± ) 
X dk^ • [b) 

From the twist expansion point of view the series Eq.fjHJ) includes the dominant part of 
all-twist T = 2 + 2n contributions ( q (2n FF) q) , (F (2n FF) F) . Explicit resummation 



technique for © was first given in [20j. The interpretation of our result is simple: it folds 
the Poisson distribution of multiple elastic scatterings with mean opacity \ = (L) /A = n, 
which depends on the quark/gluon composition of the target and the projectile, with the 
probability for transverse deflection per hit given by the normalized elastic scattering cross 
section (l/a e i) da e i/d 2 c[i±. For proton-nucleus collisions using a sharp sphere approxima- 
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tion (T A (b) = 3A^/R 2 A - b 2 /(2nR A )) and assuming equal scattering probability along the 
incident parton trajectory 

{L) = Jd 2 bT A (b) = l RA ' (7) 

where we take R A = 1.2v4 1//3 fm for the nuclear radius. One of the important conse- 
quences 0] of Eq. © is the longitudinal momentum backward shift 

-A*,, = ^ ± (8) 

that couples to the mean fcj_-broadening ( Ak"j_ ) pA = /i 2 x£- For the Gaussian approximation 
to multiple elastic scatterings, £ = const. = 0(1). Beyond this analytically convenient 
ansatz the power law tails of the distribution of the scattered partons, resulting from harder 
fluctuations along the projectile path, lead to a logarithmic enhancement of the mean k±- 
broadening. This result is the analogue of the corrections to the Moliere multiple collision 
series in QED. We use £ = ln(l + cp\) as in Ref. with c ~ 0.14 GeV~ 2 . 

To implement initial state elastic and radiative energy loss we focus on the large Q 2 ~ p\ 
partonic subprocess ab — > cd, k a , kb being the initial momenta involved in the hard part 
da ab ^ cd /dt of Eq.©. If a,b have lost fractions e a , a = a, b of their longitudinal momenta 
according to a probability distributions P a (e), at asymptotic t = — oo 

k~a = ^ ~k a , fa/p{. x a)Q ) > J de a P a (e a ) f a /p ^X a = j ~iQ J 9(%a — 1) • (9) 

Eq.® provides a simple modification to the factorized pQCD hadron production formalism. 
For bremsstrahlung processes P a (e) are sensitive to multiple gluon emission |9j. For the 
simpler case of mean energy loss P(e) = 8(e — (Ak )/k ). More specifically, for the small 
elastic longitudinal shift that we consider here 

P a (e) = s(e- , f a/p (x a , Q 2 ) - f a/p (x a + Q 2 ) 0(x a < 1) . (10) 

The observable effect of Eqs. (J9I10)) can be very different for valence quarks, sea quarks, and 
gluons due to the different x-dependence of the PDFs. 

Two fits to the yfs = 27.4, 38.8 GeV p+W/p+Be data are given in Fig.^ in the first case 
we use the EKS'98 shadowing function S a /A(x a , Q 2 ) to model the antishadowing/EMC mod- 
ification of the nuclear target. The changing quark/gluon projectile composition, reflected 
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FIG. 1: The ratio of A-scaled p + W/p + Be data on 7r + , 7r _ production at y/s = 27 A, 38.8 GeV 
from [1J|. Calculations are for |(vr + + 7r _ ): top/bottom panel include/do not include anti- 
shadowing/EMC effect as described in the text. The anticipated y/s = 200 GeV p + W/p + Be 
ratio is shown to illustrate the energy dependence of the Cronin effect. 

in its y^-dependent average color charge, enters the opacity through the elastic scattering 
cross section: 



xWs) = a elP (L) ex {C F F^, Q) + C A F^, Q)) C T 



(11) 



i.e. in a way analogous to Eq.fjHJ). In Eq. Iflljl Ct is the mean color charge of the target. The 
variation of the opacity with the center of mass energy can thus be inferred from Table H] and 
is given by x(v^) = • x(27.4GeV) for fixed Ct- While quark antishadowing/EMC 



effect are constrained by data, the modifications to the gluon PDFs are quite uncertain. In 
the second calculation we do not include nuclear shadowing. Instead, in Eq. (JTTj) we also take 
a changing average color charge for the target: Ct(>/s) = ( CpF^y/s, Q) + CAFg(y/s, Q) ). 
In this case the variation of the effective opacity with y/s is given by x(V$) = A 2 (a/s) • 
x(27.4 GeV). Both calculations include initial state mean elastic energy loss. With the 
definition of (L) given in Eq.(0) the transport coefficients of cold nuclear matter are fixed 



to fi 2 /X q w 0.06 GeV 2 /fm, /i 2 /A g = (C A /C F ) fi 2 /X q « 0.14 GeV 2 /fm. The modest ~ 20% 
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increase relative to previous estimates (/x 2 /A g ~ 0.05 GeV 2 /fm p) compensates for the 
longitudinal backward momentum shift, Eq.(jHJ), and provides comparable reproduction of 
existing low energy p + A results. This can be seen by comparing the top panel of Fig. Q to 
Fig. 1 from Ref. [5j. The bottom panel of Fig. ^ seems to show better agreement between 
data and theory in terms of improved reproduction of the Cronin peak and separation of the 
nuclear modification ratios versus y/s at high p T . We continue to explore both possibilities 
but the calculations without strong antishadowing/EMC effect are also expected to give 
better results at RHIC energies. 

The systematic reduction of the Cronin enhancement /suppression with increasing y/s has 
its natural explanation in the picture of multiple elastic scattering on a nuclear target 2(| 
(see Fig-HJ. With the opacity x changing only slightly with y/s the dominant effect becomes 
the significant hardening of the hadronic spectra as predicted by perturbative QCD. This 
reduces the observable consequences of otherwise identical transverse momentum kicks. 

We have also checked that the flavor dependence of the Cronin effect - the larger enhance- 
ment/suppression for kaons (and also protons) - cannot be reproduced within the accur acy 
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of the calculation. Studies have so far been focused on the baryon enhancement puzzle 
but the results are inconclusive. We note that improved perturbative techniques are nec- 
essary to address the Cronin effect for massive hadrons and extend the calculations to the 
Pt < 2 GeV region. 



IV. IMPACT PARAMETER AND RAPIDITY DEPENDENCE OF THE CRONIN 
EFFECT 

The y/s = 200 GeV results for the Cronin effect in d + Au reactions at midrapidity are 
shown in Fig. At y = for the corresponding 2 x 10 -2 < xt < 10 _1 values the shadowing 
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FIG. 2: The Cronin effect, represented by the nuclear modification ratio RdAu(pr), for \{tt + + vr ) 
is shown for the total invariant cross section and impact parameters 6 = 0,6 fm. 

region of S a /A(x, Q 2 ) is never reached. However, the relative contribution of antishadowing is 
seen to be large for > 3 GeV. We find that without significant x-dependent modification 
of the nuclear target (i.e. the case that gives better agreement with existing data from Fig.Q) 
the Cronin enhancement reaches 1.2 with peak position (RdAu max) in the pt — 3 — 4 GeV 
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FIG. 3: Illustration of the difference between the computed small ~ 5 — 10% Cronin enhancement 



at pt > 5 GeV in central reactions and the argued |ll( factor of 30 — 50% suppression. The 







estimated \/ Raa{pt) for \{h + + h ) and ir° from preliminary STAR and PHENIX data [lj and 
the Np ar t/2 scaling for b = 0, 6 fm are shown. 

range, a 30% drop relative to the strong gluon antishadowing scenario. In going from central 
(6 = fm) to peripheral (6 = 6 fm) reactions the magnitude of the Cronin effect is reduced 
by a factor ~ 2.5. The main uncertainty in computing RdAipr) for peripheral reactions 
comes from the large deuteron radius, leading to a distribution in impact parameter of its 
partonic constituents. The corresponding spread in opacity x = \F^a ~ b 2 /\ grows with 
increasing b and using the average value x(b) in Eqs. (j6l8|) may not be a good approximation. 

Midrapidity measurements accessible to all experiments, especially for central reactions, 
are the prime candidate that can distinguish between different hadron production scenarios. 
Glauber model calculations ^ that account for the large deuteron radius show a factor ~ 2 
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difference between jV& in and N part /2 (not N part \), the two possible scalings relative to p +p, 
over a wide range of impact parameters. Participant scaling at moderate and high px (up to 
Pt = 8 — 10 GeV as discussed in jllj]) will correspond to RdAuipr) — 0-5 as shown in Fig. 01 



The gluon saturation model suggests RdAu(pr) = \/ Ra u Au{pt) , i-e. it will in fact lead to the 
same conclusion since the suppression in Au + Au reactions is 3-5 fold. This is demonstrated 
in Fig.|3]through the estimated \/ RauAu(pt) for inclusive charged hadrons and neutral pions 
from preliminary STAR and PHENIX data [lj]. We emphasize that the px dependence of 
the nuclear modification ratio in d + Au is different as well. Final state hadron absorption 
mechanisms [l]J are also expected to result in suppression in d + Au, albeit smaller and 
difficult to quantify analytically. In contrast, the calculations in Figs. 121 and El show small 
enhancement in the specified pt range even with initial state elastic energy loss. Predictions 
that have quantitatively addressed the y = 0, pt > 2 GeV Cronin effect at RHIC from the 
standpoint of perturbative QCD and multiple elastic scatterings give similar results j^Q- 
Experimental measurements in d + Au will enable a critical test to differentiate between 
competing interpretations of Au + Au data. 

Fig. 13] shows the predicted Cronin effect at forward (y = +3) and backward (y = —3) 
rapidities near the edge of BRAHMS acceptance. In the following calculations the direction 
of the deuteron beam is chosen to be positive. At forward rapidity, which probes the small- 
x region of the nucleus, the effect of shadowing is very small. Similarly, turning off elastic 
energy loss does not significantly change the shape and magnitude of RdAuipr)- Relative 
to y = 0, both the peak and intercept (RdAu = 1) positions are shifted to slightly higher 
transverse momenta. The most distinct prediction at y — +3 is the significantly flatter 
Cronin enhancement region that extends to high p?. At forward/backward rapidities hadron 
spectra are much steeper than at y — 0. This explains the markedly broader range where 
the effect of transverse momentum kicks is observed. 

At backward y = —3 rapidity, shown at the bottom panel of Fig. HJ there is no Cronin en- 
hancement since the partons from the nucleus do not scatter multiply on the deuteron. How- 
ever, in the px range shown data probes the large x > 0.4 nuclear modification S a /A.(x, Q 2 ) 
and may provide complementary information about the EMC effect. The backward ra- 
pidity region is also more sensitive to initial state energy loss. Switching it off changes 
RdAuipr) from moderate ~ 25% suppression to an almost perfect binary collision scaling 
(RdAuipr) - 1) for p T > 3 GeV. 



11 





T3 



1.5 



1.4 



< 
< 



-< 1.3 
1.2 



- Y = 



< 

"D 

^ 1.1 



1.0 



— ' 1 ' — 

- Q.5(n+n), EKS'98 shadowing 
— 0.5(it + +Jt ), no shadowing 



1 ' 1 1 

All panels: Total invariant cross section 



2 2 2 

<k_ > = <k_ > + <Ak_ > . 

T T pp T pA 



I 





1.6 


CL 


1.5 


CL 




1 .4 








< 
< 


1.3 






< 


1.2 






< 
"D 


1.1 









1.0 




0.9 




1.2 


CL 
CL 


1.1 


D 




T3 


1.0 


< 
< 




< 


0.9 


< 


0.8 


"D 

D 




T3 


0.7 



"1 1 1 — 

0.5(it + +Jt ), EKS'98 shadowing 
0.5(jt + +jO, no shadowing 



Y = +3 



All panels: perfect binary scaling 



0.6 



i 1 1 — 

0.5(ji + +jO, EKS'98 shadowing 
0.5(it + +Jt"), no shadowing 

I ' 



. J 



T 



Upper band: no initial state E-loss 
I I i 



_L 



4 

p T [GeV] 



Y = -3 



r 



FIG. 4: Rapidity dependence of the Cronin effect in d + Au reactions at ^fs 
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200 GeV with 



and without antishadowing/EMC effect. The result of switching off elastic energy loss is also shown 
via the upper bands for y = ±3. 
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V. CONCLUSIONS 



In this letter we have derived a technique of incorporating initial state parton energy 
loss in perturbative calculations. Midrapidity moderate- and high-p^ hadron spectra show 
little or no sensitivity to this effect since the longitudinal momentum backward shifts can 
be compensated by a small increase in the estimated cold nuclear transport coefficients to 
obtain an equally good description of existing p + A data. Detecting initial state energy 
loss may be easier in the small- x region of the deuteron where the steepening of the parton 
distribution functions would tend to amplify the effect. 

In this letter we have given predictions for the centrality and rapidity dependence, see 
Figs. 121 and EJ of the moderate and high transverse momentum (p? — 2 GeV) Cronin effect 
in d + Au reactions at RHIC. They follow the systematic approach to the computation of the 
nuclear modification factors outlined in Ref. p. At midrapidity we find a small ~ 20 — 30% 
Cronin enhancement that decreases in going from central to peripheral reactions and at 
high pt- Forward and backward (y = ±3) rapidity regions are predicted to have markedly 
different Rcla{pt)- Upcoming central d + Au data will provide a decisive test to distinguish 
between suppression and enhancement dynamical models 0, [ill Q|> see Fig. El of initial 
state effects in nuclear collisions. 
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